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Equation of state for systems with particles self-assembling into aggregates is derived within a mesoscopic 
theory combining density functional and field-theoretic approaches. We focus on the effect of mesoscopic fluc- 
tuations in the disordered phase. The pressure - volume fraction isotherms are calculated explicitly for two 
forms of the short-range attraction long-range repulsion potential. Mesoscopic fluctuations lead to an increased 
pressure in each case, except for very small volume fractions. When large clusters are formed, the mechanical 
instability of the system is present at much higher temperature than found in mean-field approximation. In this 
case phase separation competes with the formation of periodic phases (colloidal crystals). In the case of small 
clusters, no mechanical instability associated with separation into dilute and dense phases appears. 
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1. Introduction 

Recent experimental theoretical liMllll and simulation @, [lil - tlill studies reveal that in many 
systems with competing interactions, clusters or aggregates of different sizes and shapes are formed. 
These objects in certain thermodynamic states can form ordered structures in space j^LEHlil]. Notable 
examples include charged globular proteins in water UHE3E3l> an d mixtures of small nonadsorbing 
polymers with charged colloids or micelles Interactions (or in fact effective interactions) in the 

latter systems can be described by the model potential consisting of short-range attraction, resulting from 
solvophobic or depletion interactions, and long-range repulsion, resulting from screened electrostatic 
potential (SALR potential). 

Systems containing clusters or aggregates are inhomogeneous on the length scale associated with the 
average size of the aggregates and average distance between them. The corresponding length scale of 
the inhomogeneities is significantly larger than the size of the particles. Fluctuations on the mesoscopic 
length scale corresponding to displacements of the aggregates have an important impact on the grand 
potential, and thus on the equation of state (EOS). Derivation of an accurate EOS for inhomogeneous sys- 
tems is less trivial than in the case of homogeneous systems, since it is necessary to perform summation 
over different spatial distributions of the clusters and over all deformations of them. 

Contribution to the grand potential associated with mesoscopic fluctuations can be calculated in the 
field-theoretic approach 0>[ll|]. In principle, this contribution can be obtained in the perturbation expan- 
sion in terms of Feynman diagrams. In practice, an approximate result can be analytically obtained in the 
self-consistent Hartree approximation A formal expression for the fluctuation contribution to 

the grand potential has been derived in references I5ll6lll9l1. However, its explicit form with the chemical 
potential expressed in terms of temperature and density has not been determined yet. The EOS isotherms 
for various forms of the SALR potential were not analyzed, and the effect of the mesoscopic fluctuations 
on pressure remains an open question. 
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It is important to note that various forms of the SALR potential are associated with different proper- 
ties of the systems. Depending on the ratios between the strengths and ranges of the attractive and re- 
pulsive parts of the potential, separation into uniform phases, formation of clusters of various sizes and 
shapes (globules, cylinders, slabs) in the so-called microsegregation, or isolated individual particles may 
occur. In certain conditions, the clusters can be periodically distributed in space in the periodic phases 
whose densities are smaller than the density of the liquid phase (i|,[lll[l2|]- Possible types of the phase 
diagram for different SALR potentials are shown in figure Q] (see also references lillillioll). Properties 



Figure 1. Types of possible phase diagrams for the SALR potentials, shown schematically. From the left 
to the right panel the role of the repulsion increases, 'p' indicates the stability region of the periodic 
phases of different symmetries. The disordered fluid close to the stability region of the periodic phases is 
inhomogeneous, i.e. clusters are formed, but their positions are correlated only at short distances. In this 
work we are interested in systems exhibiting the phase behaviour shown in the right (System 1) and the 
central (System 2) panel. 

of the disordered phase can be influenced by the periodic phases for thermodynamic states close to the 
stability of the latter. We expect that the disordered phase, although the long-range order is absent, is 
inhomogeneous on the mesoscopic length scale and resembles 'molten periodic phases'. In this respect 
the inhomogeneous disordered phase is similar to microemulsion which can be interpreted as molten 
lyotropic liquid crystal. 

In this work we focus on the stable or metastable disordered inhomogeneous phase in which clusters 
are formed, but they do not form any ordered structure. We derive the EOS with the contribution from 
mesoscopic length-scale fluctuations included. We calculate the explicit form of the EOS for two represen- 
tative examples of the SALR potential within the self-consistent Hartree approximation. The first system 
corresponds to the formation of small clusters, and the gas-liquid separation is unstable for all tempera- 
tures in the mean-field (MF) approximation (figure[l] right panel). In the second system, large clusters are 
formed. The gas-liquid separation is present in this system as a stable or a metastable transition for low 
temperatures (figure [1] central panel). We shall compare the effect of mesoscopic fluctuations in these 
two cases on the isotherms P(0» where ( is the volume fraction of particles and P is pressure. 

In the next section we briefly summarize the mesoscopic approach. In section 3 the EOS is obtained 
by two methods. In section 3.1 we consider mesoscopic fluctuations about the average volume fraction, 
while in section 3.2 fluctuations about the most probable volume fraction are included. The formulas de- 
rived in section 3 are evaluated for the two versions of the SALR potential in section 4. We obtain a com- 
pletely different effect of the mesoscopic fluctuations in these two cases. The two approaches (sections 3.1 
and 3.2) yield very close results provided that the relative fluctuation contribution to the average volume 
fraction is small. For larger fluctuation-induced shifts of the volume fraction, only qualitative agreement 
of the two methods is obtained. Short summary is presented in section 5. 



We consider a local volume fraction of particles, i.e. the microscopic volume fraction averaged over 
mesoscopic regions, as an order parameter |6]. The corresponding mesoscopic volume fraction varies 
on a length scale larger than the size of the particles, and the characteristic size of inhomogeneities is 
the upper limit for the mesoscopic length scale. A particular form of the mesoscopic volume fraction 
can be considered as a constraint on the microscopic states. The corresponding mesostate is a subset of 
microstates compatible with the imposed constraint. The mesoscopic volume fraction (or the mesostate) 



T 



T 



T 




2. Short summary of the mesoscopic description 
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was defined in references 00]. For a one-component case we fix the mesoscopic length scale R 3= a/2 
and consider spheres Sr (r) of radius R and centers at r that cover the whole volume V of the system. We 
define the mesoscopic volume fraction at r by 

Cto:=jjr / C(r',^), (2.1) 
r'ES fi (r) 

where V$ = 4tiR 3 /3, and the microscopic volume fraction in the microstate Jt = {{r ! } i= i j jvl is defined 

by 

N to \ 

Ur,J£):=Y,0 (— -|r-r f | ), (2.2) 

where Q{r) is the Heaviside unit step function. The microscopic volume fraction is equal to 1 at points 
that are inside one of the hard spheres, and zero otherwise. Integrated over the system volume, it yields 
the volume occupied by the particles. The mesoscopic volume fraction at r is equal to the fraction of 
the volume of the sphere Sr(t) that is occupied by the particles. Note that f(r) takes the same value if 
one particle is entirely included in this sphere, independently of the precise position of its centre. Thus, 
C(r) gives less precise information on the distribution of particles than £(r). In the disordered phase f(r) 
is independent of r and equals the fraction of the total volume that is occupied by the particles. The 
mesostate can be imagined as a fixed distribution of centers of clusters, with arbitrary distribution of 
particles within the clusters, and small modifications of their shapes. Probability of the mesostate £ is 
given by 

e -/JGcoK] 

pl(] = = , (2.3) 



where 

w 

D( e -PQcol(]_ (2.4) 



The functional integral / D( in <2.4) is over all mesostates, 

Q co = U[0-TS[(]-nN[(], (2.5) 

where U, S, N are the internal energy, entropy and the number of molecules respectively in the system 
with the constraint of compatibility with the mesostate [(] imposed on the microscopic volume fractions. 
U is given by the well known expression 

U[0 = \J J v co(n ~ r 2 )C(ri)C(r 2 ), (2.6) 

ri r 2 

where for spherically symmetric interactions 

^co(ri-r 2 ) = ^(ri 2 )gco(ri-r 2 ), (2.7) 

r\2- lri-r 2 |, v 2 V{r\2) is the interaction potential, v — na 3 !6 is the volume of the particle, andg co (ri-r 2 ) 
is the microscopic pair correlation function for the microscopic volume fraction in the system with the 
constraint of compatibility with the mesostate imposed on the microscopic states. The grand potential 
can be written in the form 

)3nK] log H fluc , (2.8) 



where 



E Bac = flkf>e~ pHltai ^, (2.9) 



y [ f^l 



H nilc [i:,<p] = n co i(+(f)]-n co [(]= V I ... I ^-(p(r 1 )...(f>(T n ). (2.io) 
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The average mesoscopic volume fraction, (, corresponds to the minimum of Q, and must satisfy the equa- 
tion 

fiffWfi , / g(/?Hflac) \_ (2n) 



where the averaging is over the fields cp with the probability oc exp{-/3Hfi uc [(,(p]). Note that when 
Cf° [(] = for odd n, then the second term on the LHS in 12.111 vanishes, and the average volume fraction 
coincides with the most probable volume fraction (q given by 



5CW 



= 0. (2.12) 



By contrast, when C™[f] # for odd then C ^ Co- 
in order to evaluate the fluctuation contribution to Q.[(] we decompose Hr uc [( , cp] into two parts 

H nuc [(,c/)] = J? G [(,4>] + A,je>[(,4>], (2.13) 

where in the disordered phase 

•afblf.fl = ^ J 0(k)C 2 (k,O0(-k), (2.14) 

k 

and Q (k, () is the Fourier transform of 

c 2 (ri 2 ,fl = — H . '7 (2.15) 

5<(ri)5f(r 2 ) 

The above function calculated for £ = Cis related to the direct correlation function j2oll . 
Assuming AJ£ « J£q, we obtain la Hill 

pniQapOcolft-logJ D<pe-^ + {pAJe) G +0{{pAJe) 2 G ), (2.16) 

where (. . .)g denotes the averaging with the Gaussian Boltzmann factor oc e~^ G . 

Since P--Q./V, approximate EOS can be obtained from d2.16> calculated for ( satisfying I2.11K when 
^co 10 (see H2.5t ) is known. The chemical potential in H2.16t should be expressed in terms of T and (; its 
form as a function of T and ( can be determined from equation 12.11) . In order to evaluate the second 
term in equation <2.11> . we need approximate forms of the correlation functions. In the lowest order 
approximation, it is necessary to determine 

<<p(k)0(-k)> = G 2 (k,0 = l/C 2 {k t (). (2.17) 



3. Approximate results for the fluctuation contributions to the EOS, 
density and chemical potential 

In this section we derive an explicit form of the EOS P{(, T), the density shift, A( = (-(0, and the chem- 
ical potential T) under the following assumptions: (i) local density approximation for S[(] and (ii) the 
lowest-order approximation for the second term in <2.11> . In the local density approximation we have 



TSIO 



h(CW), 



(3.1) 



where fh(0 is the free-energy density of the hard-sphere system with dimensionless density p* 
We assume the Percus-Yevick approximation 



PM0 = P*m P *)-p*+p* 



3C(2-Q 
2(1 -0 2 



ln(l-C) 



6(1 71. 



(3.2) 
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In the local density approximation C{j° [(] are just functions of ( in the disordered phase, and we can 
simplify the notation, introducing 

For n > 2 we have 

C™[0=A n iO, (3.4) 

whereas for n = 2 

C?ik,0 = PV co {k) + A 2 {0, (3.5) 
where in the disordered phase C^tfc.O is the Fourier transform of the function 

i!£ML. (3 . 6 ) 
2 5f(n)5C(r 2 ) 



Equations 13.2) and (276) define the functional Q co for a given form of V co (equation (277)). 

From equation (2773) it follows that the most probable fluctuations correspond to the wavenumbers 
k = fc D for which V co {k) assumes the minimum, and the inhomogeneities on the length scale Inlk^ are 
energetically favored when Vco(^b) < 0- In this work we focus on the effect of the self-assembly into 
aggregates. Therefore, we restrict our attention to Vco (A:) which assumes the minimum for k b > 0, and 
Vco(^b) < SEHl- Since the fluctuations with the wavenumber k a k b are most probable, they yield 
the main fluctuation contribution to the grand potential 12.16) . For such fluctuations we can make the 
approximation 

pV C0 {k) « pv C0 {k h ) + pv£Hk h ){k- fc b ) 2 /2 + ... (3.7) 

As the energy scale we choose the excess energy associated with the fluctuations having unit amplitude 
and the wavenumber k^, and introduce the notation 

P* = UT* =j3|Vco(fcb)l, (3.8) 



* = _ V ™ {kb) (3.9) 



2|V C0 (fc b )| 
and 

V ; = S^L. (3.io) 

3.1. Fluctuations around the average volume fraction 

In this subsection we consider fluctuations about the average value (. We shall first determine the 
chemical potential as a function of ( and T* from I2.ll> . In order to calculate the second term in I I2.11) , 
we assume that relevant fluctuations are of small amplitudes, and truncate the expansion in 12.10) at the 
fourth-order term. Next we insert the derivative with respect to ( of the RHS of equation (2.10) truncated 
at the quadratic term in <p, and we obtain from 12.11) and 12.5) an approximate equation for the rescaled 
chemical potential, p, = 6[il(n<j 3 ), of the form 

/5/2*/3A MF (0 + : ^»(0, (3.1D 

where 

Pi* MF (0=pv co m+pf{ i (0 p.12) 
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and the last term in (3.11) is the fluctuation contribution with 

<S(.0 = jG2Vc,0 = G2(Q,Q. (3.13) 

k 

The same expression can be obtained from <(/>} = 0, when <2.13t . 1 12. 14) and the approximation 

e -/Sfffiuc _ e -PH G 1 1 _ p AH + [ AH 2j 1 (3 14) 

are used. The approximation <3.11> is valid as long as the correction term is not larger than the MF result. 
When considering particular cases we shall verify if this is the case. The fluctuation contribution to the 
direct correlation function 42.15) is obtained by calculating the second derivative of the second term on 
the RHS of <2.8) with respect to (. In the consistent approximation we insert in the obtained expression 



the appropriate derivatives of H^ uc ^equation 12.10) ) with the expansion in (p truncated at the second 
order. The result is given by 



of ti« m . leg 

(III EI EI 



C 2 (fc,0«C 2 co (fc,0 + ^^(0. (3.15) 

Equations 13.15) , 12.171 and 13.131 should be solved self-consistently. 

The fluctuation induced shift of the volume fraction, A( - ( - (o, can be obtained from H2.ll) by ex- 
panding the first term on the LHS about Co. 

C T (0 ,Co)AC + £ = _/*a^\. (3.16) 

For small A( we can truncate the expansion in J3.16) at the first term. When the RHS in 13.161 is approxi- 
mated as in the calculation of ft from J2.ll) , we obtain the result 

a _ ^(Co) ^ a __ A3K0) 

2C 2 «Uo) 2C 2 {0,(o) + A 3 (( )^(( ) 

We used the approximations: A 3 (Q = A 3 {( ) + A 4 (( )AC + 0(AC 2 ), #(0 = <^(Co) + ^'{()A( + 0{A( 2 ) and 
equation 13.15) . 

For the potential given in @77) the approximate form of ^ is (ijlliHililiill 

2a\/f* 

^0 = -^-. (3.18) 



where 



a= ^= (3.19) 



and Z(0 = \fC 2 {k h ,Q. The above approximation is valid for C 2 (k b ,0 « P*v*kl d [H, [H 0]. The 
equation 13.15) for k=k\ } takes the form 

Z(0 3 = Z(0££°(*b,0 + ^(Oflv^, (3.20) 

and the explicit expression for Z is 

Z(0 = — — + — ? — (3.21) 

6 W(0 

with 

W(0 = j 108^4(0 a \/f* + 12 \J-l2Cf (fc b , 3 + 81 [ ^ CO a vT*] 2 1 . (3.22) 
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The fluctuation contribution in equation 12.16) for the approximations (3.15M3T201 was calculated in 
references |5j|l9|,|2l|,|22|], and has the form 



pn [f] « ^n co iQ + 2a VWz{0 v - \ v. 

o 

Taking into account (3. Ill and (3?2) , we obtain from 13.23) the explicit form of the EOS 

PP(0 = PP MF {0 + F{(, T*), 



where 



and 



2 d( 

aVT*A 3 tO( 



(wf 



F{(,T*) 



2aVT*Z(0 + 



2A 3 {Q(-q o {k h ,Q 
Z(0 



a A 4 (()T* 
2Z(0 2 



•3Z(0 



(3.23) 



(3.24) 



(3.25) 



(3.26) 



The second equality in 13.25) is valid for the PY approximation for f^. In order to obtain the last equality 
in <3.26t . equation 13.20) was used. 



3.2. Fluctuations around the most probable volume fraction 

In the previous subsection we considered fluctuations about the average value, which in general dif- 
fers from the most probable value of the volume fraction. In principle, it is possible to consider equations 
analogous to 12.8) and 12.10) . but with ( replaced by (q. In this new approach equation 12.12) is satisfied, 
and thus the expansion in <2A0l starts with n = 2 (CJ° l( ] = 0). On the other hand, (</>>= C - Co + 0. The 
results obtained in the two approaches — with included fluctuations around the average value or around 
the most probable value — should be the same in the exact theory. However, when the fluctuation contri- 
bution is obtained in an approximate theory, the results may depend on the validity of the assumptions 
made in the two approaches. In this section we derive an alternative version of the EOS, based on the con- 
tribution from the fluctuations around the most probable value. From 12.12) we obtain for the chemical 
potential 

A«o) = /i MF «o) (3.27) 

with fi MF given in equation 13.12) . We consider 12.13) with ( replaced by (o, and the approximation 13.14) . 
In the above, Hq is given in equation 12.14) with Qtfc.Co) = \IG2(k,(o), where G2(fc,<To) = (0(k)0(-k)), 
and 

AH = \ff ^(r')[cr(r'-r",( o )-C 2 (r'-r",Co)]0(r") + / ^fV) 3 + f M^0(r) 4 + ... 
r' r" r r 

(3.28) 

Note that the function G2(fc,<To) defined here differs from the correlation function, because in this case 
<0(r)> + 0. Taking into account that for H G of the form <Q3) there holds / D(pe~P HG (p 2n+1 = 0, we obtain 
the expression 



Finally, the lowest-order result is 



I = ft, + <0) » ft - n ^f n °l . g(Co). (3.30) 
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At the same level of approximation C2(fc,Co) is given in equation 13.15) , except that all quantities are 
calculated at Co which satisfies 12.121 rather than I2.11> . This can be verified by a direct calculation 
of <(/>(ri)(£(r2)} with the help of 13.14) and 13.28) , in an approximation analogous to <3.29> (see refer- 
ence HI). 

The above shift of the volume fraction differs from I3.17I . because instead of ^(C), there appears 
^(Co)- The dependence of p. on the average volume fraction is given in equations J3.27I and 13.30) . with 
eliminated Co- 
in order to evaluate the EOS, we consider an equation analogous to 1 13. 23t . with fi co [Co] calculated at 
its minimum ( = Co- The EOS takes the form 

j3P(Co) = ,0P MF (Co) + Fo(Co. r*), (3.31) 

where fiP M¥ is defined in (3251 , Co satisfies fTTZi , and 

— a 2 A 4 (.( )T* 
F {(o,T*) = -2av r T T Z(( ) + * ' ' . (3.32) 

2Z(Co) 

The dependence of P on the average volume fraction ( is given by parametric equations 13.311 with <3.32> 
and fOuV 

The approximate theory developed in this section is valid for small AC since we assumed that the 
relevant fluctuations are small and truncated the expansion in <2.10> at the term oc cp 4 . Moreover, to 
evaluate AC we neglected the terms of the order 0(AC 2 ). We may expect that if we obtain large AC and 
large discrepancies between the results obtained by the two methods, then the approximate theory is not 
sufficiently accurate. 



3.3. Comparison between the two methods 

Let us focus on the chemical potential, and compare the two expressions, equations 13.11) and 13.27) 
where ( and Co satisfy equations J2.ll> and 12.12) . respectively. We expand the RHS in equation 13.11) 
about Co. 

fi a j8V- co (0)Co + A/h(Co) + AC 

+ ^Wo) + 0(AC 2 ). (3.33) 

From 13.171 , 13. 15) and 13. 271 we obtain an equality of the two expressions for the chemical potential to 
the linear order in AC, when AC is given in 13.171 . 

Similarly, to compare the two expressions for the EOS, equations 13.24) and 13.31) , we expand the 
RHS of equation 13.24) about Co to the linear order in AC- Taking into account 13.17) and 13.20) , we arrive 
at equation 13.31) , up to the terms proportional to As. The latter are disregarded in an approximation 
consistent with the c/> 4 theory for the fluctuation contribution considered in this work. For relatively 
large AC, when the terms beyond the linear order become important, discrepancies between the results 
obtained by the two methods should be expected. 



^4 (Co) 

jSKo(O) + A 2 (Co) + -y^CCo) 



4. Explicit results for two model potentials 

In this subsection we shall compare the expressions for the chemical potential and for the pressure 
obtained by the two approaches for two systems showing a qualitatively different behavior. We shall 
evaluate the EOS 13.24) for the representative model SALR potential, 



V co ir): 



-z 2 r 



6{r-l), 



(4.1) 
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where z\ is the inverse range in a 1 units. The function 6{r - 1) is a very crude approximation for the 
pair distribution function. In Fourier representation, the above SALR potential takes the form 



Vco(fc) = 47T 



srfze 7,2 I sinfc 

4^{ Z2 — +cosk 



zj + k 2 



Zl I sink 

Z\ — h COSfc 



(4.2) 



We choose two sets of parameters, considered in reference fllll in the context of most probable inho- 
mogeneous structures 



System 1 : 
System 2 : 



1, sl 2 = 0.05, zi = 3, z 2 = 0.5; 
1, sl z = 0.2, zi = 1, z 2 = 0.5. 



The relevant parameters, fcb, and Vg, (see d3.9t and <3.10) ) take the values: v\ ~ 3.02 and 



System 1 : 
System 2 : 



fc b a 1.79, 
k h a 0.6089, 



Vq a -30.145; 
Vq a 0.035. 



(4.3) 



(4.4) 



The two potentials in Fourier representation are shown in figure [2] In the first system small clusters 
are formed, since 2x1^ is small. Moreover, V co {0) > 0, and the clusters repel each other. The gas-liquid 
separation is entirely suppressed due to the very short range of the attractive part of the potential. In 
the second system large clusters are formed, and V co (0) < (the clusters attract each other). Therefore, 
in MF, the metastable separation into disordered low- and high density phases occurs at low tempera- 
ture. Simulation results show that when large clusters are formed, gas-liquid separation occurs for low 
temperatures, and periodic phases are stable at higher temperatures (9J]. 



v„(k) i 




4 6 8 10 

k 



T* 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 



to 



Figure 2. The potential V co {k) for System 1 (solid 
line) and for System 2 (dashed line). V co {k) is in 
dimensionless units, k is in units tr -1 , where a is 
the particle diameter. 



Figure 3. The universal structural line (solid) in re- 
duced units (see (3.8) ) and the metastable spinodal 
line of the separation into dilute and dense phases 
(dashed) for System 2. Note that the temperature 
scale is different from the corresponding scale in 
the theory based on p* I'll], and the scaling factor 

iS (6/7T) 2 . 



We are interested mainly in the part of the phase diagram where the homogeneous structure is less 
probable than periodic distribution of particles in space, i.e. when O c0 does not assume a minimum for 
(o = const (see J2.3) ). We stress that the most probable structure differs from the average structure due 
to mesoscopic fluctuations. Cluster formation is associated with the excess volume fraction followed by 
a depleted volume fraction in mesoscopic regions, and the most probable mesoscopic structure associ- 
ated with cluster formation is periodic. Displacements of the clusters (i.e., mesoscopic fluctuations) can 
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destroy the long-range order, though. Indeed, when temperature is sufficiently high, the average volume 
fraction ((r) takes the constant value as a result of the averaging over cluster displacements, and the dis- 
ordered inhomogeneous structure with short-range correlations of the cluster positions is found I22II23I1 . 

1 1 1 — 1 

On the other hand, for low temperatures, the ordered periodic structures are stable 1 22, 23]. The phase- 
space region where the inhomogeneous phases (with either short- or long-range order) are stable is en- 
closed by the structural line fl|2J] given by T* = 1/AziO an d shown in figure[3] The structural line is also 
referred to as A-line in literature la liol 1241291 . Note that in the reduced units (see J3.8) ) the structural 
line is universal. 




Figure 4. A^(Co) f° r System 1. Dotted, solid, 
dash-dotted, and dash lines correspond to T* = 
0.02,0.01,0.007 and 0.0008, respectively, equa- 
tions (3.30) and 13.16) are not distinguishable on 
the plot. 



Figure 5. Af (f 0) for System 2 from equation 13.301 . 
Dotted, solid, dash-dotted and dash lines cor- 
respond to T* = 0.02,0.01,0.007,0.0008, respec- 
tively. 



We first compare the change of the average 
volume fraction induced by mesoscopic fluctu- 
ations. The shift A( calculated from 1 13.30) and 
( 13.16) in System 1 is shown in figure [4] The shift 
is small for a relevant range of temperatures, and 
both formulas yield practically the same result — 
they are indistinguishable on the plot. The shift in- 
creases for a decreasing temperature. In System 2, 
the fluctuation contribution to the volume frac- 
tion is much larger than in System 1 (figure [5). As 
expected, when A£ is not very small, &(/(o ~ 0.2, 
then the two approaches yield somewhat different 
results, as shown in figure [6]for System 2. 

In the next step, we study the chemical poten- 
tial. The fluctuation contribution in System 1 is 
small, except at very small volume fractions (fig- 
ure [7}, whereas in System 2 it is substantial, and 
increases for decreasing temperature, as shown 
in figures [8] and [9] The two approaches yield sim- 
ilar results for small AC, whereas when A(/(o > 
0.25, significant discrepancy between the two ap- 
proaches is obtained. We can conclude that on the 
quantitative level the approximate theory is over- 




Figure 6. Comparison of the two expressions, equa- 
tions COol and (3361 , for A£(Co) for System 2. The 
lines with smaller \A(\ correspond to T* = 0.02; the 

, and the dotted 



solid line represents equation 13.16) 
line represents equation (3.30) . The remaining lines 
represent equation 13.16) (solid) and equation <3.30) 
(dashed) for T* =0.01. 
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Figure 7. Left panel: the relative difference, (/i-/i MF )//j MF , between fi calculated from equation 13.11) . 
and the MF approximation /i MF (3A2\ for System 1 at T* = 0.007 and ( s 0.01. Right panel: p, - £ MF for 
System 1 (solid line) and System 2 (dashed line) for 0.0005 « ( < 0.1. 



simplified for the range of T and ( for which there are significant discrepancies between the two ap- 
proaches. 

Note that since p. - ji MF is large for volume fractions ( < 0.03, for very small volume fractions our 
results are oversimplified. 

Finally, we present the isotherms obtained from d3.24t and (3.311 for the two systems in figures lToHTH 
In System 1 the pressure is much higher than found in MF, and for all temperatures it monotonously 
increases with (, as shown in figure ITOl The increased pressure associated with mesoscopic fluctuations 
may result from the repulsion between the clusters, because in this case Vco(0) > 0. 



30- 



■ v.- 



0.2 



0.1 0.2 0.3 0.4 0.5 



Figure 8. Chemical potential in MF (dotted line), 
and with the fluctuation contribution included 
according to equation (3.11) (solid) and equa- 
tion CQ7) with rOol (dashed line) for T* = 0.015 
in System 2. The volume fraction is dimensionless, 
and the chemical potential is in jikT/6 units. 



Figure 9. Chemical potential in MF (dotted line), 
and with the fluctuation contribution included 
according to equation (3.11) (solid line), equa- 
tion (3Tf\ with equation (330\ (dashed line), 
and equation (3.27) with equation (3.16) (dashed- 
dotted line) for T* = 0.007 in System 2. The vol- 
ume fraction is dimensionless, and the chemical 
potential is in nkTIG units. 
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Figure 10. The EOS isotherms for System 1 
for T* = 0.02. Dashed line is the MF re- 
sult (equation 13.25) ) and solid line represent 
equation 13.24) . and equation (3.31) with equa- 
tion 13.30) . indistinguishable on the plot. 



Figure 11. The EOS isotherms <3T24l (solid line) 
and 13.31) with equation 13.30) (dash-dotted line), 
and the MF approximation 13.25) (dashed line) for 
System 2 for T* = 0.015. 



Since in System 2 Vc O (0) < 0, a mechanical instability develops at the MF spinodal line, with the 
metastable MF critical point r c * a 0.0009. What is really interesting is that such instability appears at 
much higher temperatures due to mesoscopic fluctuations. This is in strong contrast to the fluids with 
purely attractive interactions, where density fluctuations decrease the critical temperature with respect 
to the mean-field estimate. The present case with dominant fluctuations associated with mesoscopic 
wavelengths bears some resemblance to the restricted primitive model (RPM) of ionic systems. There 
is no gas-liquid instability in the RPM at the MF level of a mesoscopic theory analogous to the one con- 
sidered here, but when the short-wavelength charge-density fluctuations are included, such instability 
appears ]24j]. One could imagine that the mesoscopic fluctuations, i.e., displacements of the clusters from 
their most probable locations lead to their coalescence when Vc O (0) < 0, and thus support the phase sep- 
aration. Note that properties of System 2 are completely different from the previously studied System 1. 




Figure 12. The EOS isotherms <3T24l (solid line) and EOT) with equation (330\ (dash-dotted line), and the 
MF approximation 13.25) (dashed line) for System 2 for T* = 0.0093. Left and right panels show large and 
small range of volume fractions, respectively. 
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Figure 13. The EOS isotherms <T24l (solid line) 
and (3.31) with equation 13.30) (dash-dotted line), 
and the MF approximation (3.251 (dashed line) for 
System 2 for T* = 0.007. 



Figure 14. The EOS isotherms f3T24> (solid line) 
and (3.31) with equation (3.30) (dash-dotted line), 
and the MF approximation (3.251 (dashed line) for 
System 2 for T* = 0.005. 



In this region of the phase diagram, the shift of the volume fraction is large. Therefore, on the quanti- 
tative level, the results are not sufficiently accurate. The inflection point on the P(() isotherm appears at 
T c * a 0.0093 or T c * = 0.007 according to ( EOT) or (3T24) with (3301 , respectively. Both temperatures, how- 
ever, are much higher than in MF. Further studies are required to verify if the separation into disordered 
inhomogeneous phases, or periodic ordering of clusters occurs. If the ordered phases are formed, the still 
open question is for which part of the phase diagram such phases are globally stable. 

Finally, let us focus on the pressure for very small volume fractions. In the fluctuation correction 
to pressure (equation d3.26» the first term comes from the fluctuation contribution to chemical poten- 
tial (see J3.ll) ). As shown in figure [7] (right panel), for very low volume fractions, our approximation is 
oversimplified, so the negative pressure is an artifact. For very low volume fractions, we should expect a 
perfect gas behavior, except that some fraction of particles should form clusters. Pressure should be pro- 
portional to the sum of the number densities of monomers and clusters. Since the number of clusters is 
smaller than the number of particles forming them, pressure should be smaller than in the corresponding 
perfect gas of isolated particles. Our theory agrees with this expectation (see figure 1121 right panel). 



5. Summary 

In this work, the effects of mesoscopic fluctuations on the average volume fraction, chemical potential 
and pressure as functions of temperature and the average volume fraction were considered within the 
framework of the mesoscopic theory | 5, 6]. We restricted our attention to a stable or metastable disor- 
dered phase. The fluctuation contribution to the quantities mentioned above was calculated in two ways. 
First, we considered the fluctuations about the average volume fraction, and derived equations 13.16) , 
1 13.11) and J3.24) for the volume fraction, chemical potential and pressure, respectively. In the second 
version, we considered the fluctuations about the most probable volume fraction, and obtained equa- 
tions (330), (321) and (331) . with ( satisfying (2A2) . The chemical potential and the EOS as functions of 
the average volume fraction are given by parametric equations 1 13.30) and d3.27) , and d3.31) , respectively. 
Our expressions are derived under the assumption that the dominant fluctuations are of small ampli- 
tudes. Consistent with the above assumption, the two methods yield the same result to a linear order in 
the fluctuation contribution to the volume fraction A( . 

The fluctuation contributions to all three quantities were explicitly calculated for two versions of the 
SALR potential. In System 1 the zeroth moment of the effective interactions is positive and small clusters 
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are formed. In System 2 the zeroth moment of the effective interactions is negative, and the clusters 
are large. We obtain nearly the same results independently of the method used when the fluctuation- 
induced shift of the volume fraction is very small (System 1). When A(/(o ~ 0.2 (low temperature in 
System 2), significant discrepancies between the two methods appear for some part of the phase diagram. 
The largest discrepancies are present when both methods yield the results that strongly deviate from the 
MF predictions. The larger is the probability of finding inhomogeneous mesoscopic states compared to 
the homogeneous distribution of particles, the stronger are the discrepancies between the two methods. 
It is in this part of the phase diagram that the periodic order may appear. We conclude that the first 
method is superior to the second one because it is easier to implement. Exact results would be necessary 
to get a comparison between the accuracy achieved by these methods. 

We have found that mesoscopic fluctuations play a very important role and lead to a significant 
change of the chemical potential and pressure. The larger is the probability of finding the inhomo- 
geneities, i.e., the further away from the structural line on the low-T side of it (figure [2), the larger is 
the role of fluctuations. When small clusters are formed (System 1 in section 4), the fluctuation contri- 
bution to pressure increases monotonously with an increasing volume fraction. By contrast, for large 
clusters (System 2 in section 4) the fluctuation contribution to pressure is nonmonotonous; it is negligi- 
ble for small as well as for large volume fractions, whereas for intermediate volume fractions it is large 
and increases with a decreasing temperature (figures IT0HT31 . Moreover, an inflection point at the pres- 
sure - volume fraction isotherm appears at the temperature and volume fraction both much larger than 
found in MF. Further studies are required to verify if the periodically ordered cluster phases are stable 
in System 2, or phase separation occurs due to the mechanical instability. Possible scenarios are: (i) phase 
separation at low T, and periodic structures at higher T, or (ii) the phase separation is only metastable, 
and finally (iii) the periodically ordered phases are only metastable. For System 1, the phase separation 
is not expected. 

In addition to the assumptions discussed earlier, we make an approximation concerning the form 
of the effective potential V co . Note that in equation <4.1> we assumed that the pair distribution function 
vanishes for r < 1 in a units. For the volume fraction, this is a poor approximation, and quantitative 
results for the structural line depend on the form of the pair distribution function (regularization of the 
potential jUl). 

In the future studies, the EOS for periodically ordered cluster phases should be determined in order 
to find the phase diagram. Our results indicate that despite the universal properties of the dependence 
of the most probable structures on the thermodynamic state, the effect of fluctuations on the average 
distribution of particles may depend on the form of the interaction potential, especially on the sign of the 
zeroth moment of the effective interactions. 
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Bn/iMB Me30ci<oniMHMX <j>/iyi<Tyai4iM Ha piBHAHHfl CTaHy 

K7iaCTep0yTB0pK)Ba/1bHI/IX CMCTeM 

A. Lrf o. nauarahP 

iHCTHTyT c|)i3MHH0T xiMii, no/ibCbKa aKafleMifl HayK, 01-224 BapiuaBa, no/ibina 
IhtCTHTyT cf>i3i/iKi<i KOHfleHCOBaHwx cudeM HaijiOHa/ibHoT aKafleMiT HayK YKpaiHn, 
syn. CBeHL(ii4bKoro, 1, 79011 JlbBiB, YKpaiHa 

PiBHAHHfl CTaHy fl/ia cudeM nacTHHOK, 140 caMOCKynnyioTbCfl b arpera™, e OTpuMaHe b paMKax Me30CK0ni- 
4hoT TeopiT, mo noeflHye MeTOfl cfiyHKLjiOHa/iy rycTWHW i TeopeTHKO-noyibOBHfi niflxifl. Mn floaiiflxyeMO Bn/inB 
Me30CKoni4Hnx cfi/iyKTyaijM y HeBnopaflKOBaHM <f)a3i. Abho o64wc/ieHO i30TepMW 'tuck - o6'eMHa nacTKa' fl/ia 
flBOx Ha6opiB napaMeipiB noTeHL(ia.ny KopoTKOCflXHe npuraraHHH n/iioc ,qa/ieK0CflXHe BiflLUTOBxyBaHHH. B ko- 
XHOMy BwnaflKy BpaxyBaHHfl Me30CK0ni4Hwx (f>/iyKTyaL(ii7i npwBOflHTb flo niflBumeHHH TUCKy, 3a bwhhtkom flyxe 
Ma/iwx o6'£mhwx 4acT0K. Ko/iw yTBoproiOTbra Be/iwKi K/iacrrepw, MexaHHHa HecTitiKicTb cucTeMw npucyTHa npn 
Ha6araT0 bwlhhx TeMnepaTypax, Hix u,e 6y/io OTpwMaHO b Ha6/iwxeHHi cepeflHboro no/ia. B i^bowiy BwnaflKy 
4>a30Be BiflOKpeM^eHHfl KOHKypye i3 c|>opMyBaHHflM nepioflw^Hnx c(>a3 (KO/ioTflHwx KpucTa/iiB). y BunaflKy Ma/iwx 
KJiacTepiB MexaHi^Ha HecrifiKiCTb, noB'fl3aHa 3 BiflOKpeM/ieHHflM b po3piflxeHy i rycry 4>a3n, He BUHUKac. 

ICmowoBi c/iOBa: Knacrepi/i, caMOCKyrmeHHH, pibhrhhr ciany, MG30CK0ntHHi <pnyKTyau,i'i 
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